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We present analytical methods (i) for calculating microwave radio re- 
fraction for negative and positive initial ray angles accounting for station 
height and (ii) for determining refraction-corrected ranges of antenna point- 
ing azimuth within which mutual interference with geostationary satellites 
in shared frequency bands is likely. 

I. INTRODUCTION 

When radio-relay and communication satellite systems share fre- 
quency bands, as they do at 4 and 6 GHz, it is necessary to impose 
restrictions on both systems so that interference is not excessive. The 
CCIR (International Radio Consultative Committee) recommends 
that radio-relay antennas maintain a specified angular separation with 
respect to the geostationary (stationary equatorial) orbit or, where 
this is not practicable, the application of power limitations to terrestrial 
radio transmitters involving reception at the satellite. While the above 
restrictions protect satellites, designers of terrestrial systems should 
be aware of possible interference into radio-relay systems from satellite 
radiation arriving at low elevation angles and close to the on-beam 
directions of receiving antennas. Because the dielectric constant of the 
earth's atmosphere varies with altitude, the radio-relay beam is not 
straight, and atmospheric refraction must be considered when computing 
the directions of radio beams for which the restrictions apply. 

l.i Simplified Exposure Model 

Figure 1 introduces the geometry of the problem and illustrates 
significant trends and limits. Radio-relay site P located at North Lati- 
tude <p degrees is shown as viewed from above the earth. An arc of the 
geostationary satellite orbit is also shown. The orbit longitude of point 
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ANTENNA POINTING 3389 

S p and the earth longitude of site P coincide numerically (that is, S p 
is in the direction of true south, PS P , observed from P). 

Two radio rays from P intercept the orbit symmetrically at points S. 
The controlling geometric relationships are based upon triangles formed 
by points S, P, and geocenter 0. Angle O-P-S is determined by the 
elevation angle of the radio-relay antenna including ray bending due 
to atmospheric refraction. A special case is depicted in Fig. 1, wherein 
unbent intercepting rays PS and also ray PS P are assumed to he in the 
local horizontal plane at P and remain tangent to the earth sphere at P. 
Thus, angles S-P-S p are always antenna pointing angles to orbit inter- 
ception referred to "true south." 

It is instructive to visualize the relationship between latitude <p at P 
and the location of orbit intercepts S, for a given fixed triangle OPS. 
As points S approach S p , the constraints imposed above require that 
the station latitude y approach a maximum latitude "visible" to the 
orbit. The resulting single pointing direction to orbit intercept is due 
south (from P to S p ). Conversely, the maximum separation between 
points S and S p obtains when $ is zero (for site P located on the equator). 
Since both intercepting rays are tangent to the equator at P, the limiting 
pointing directions are due west and due east. 

Note that a rotation in azimuth of the antenna (about the local 
vertical axis at P) between known orbit-intercept directions results in 
radio rays PS which fall below the orbit as viewed from P; rotations 
beyond these "critical azimuths" result in rays above the orbit. 

1.2 Computations 

Given the latitude and elevation angle of a microwave radio-relay 
antenna, one can calculate the pointing azimuth for which, neglecting 
refraction, the main beam axis intercepts the geostationary orbit. This 
calculation is repeated to produce screening charts like Fig. 2. Such 
charts are adequate for quickly determining a hazard condition, but 
often the true critical azimuths must be approached closely while main- 
taining tolerable interference levels. 

A graphical procedure proposed for the convenience of system plan- 
ners provides pointing angle estimates for most stations when caution 
is exercised in those steps accounting for atmospheric refraction. 1 An 
analytical technique adaptable to machine calculation is also required 
for rapid, accurate screening of large numbers of existing and proposed 
radio-relay sites for potential interference exposures. Precise evaluations 
are required for cases of unavoidable exposure. 

The method described in following sections can be used by the system 
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Fig. 2 — Screening chart, neglecting refraction. 

planner to determine critical azimuths and the ranges of azimuths to be 
avoided (accounting for atmospheric refraction) for prescribed minimum 
angular separations between main beam axes and the geostationary 
orbit. 

The Central Radio Propagation Laboratory (CRPL) Exponential 
Reference Atmosphere is adopted for the generation of microwave 
radio refraction curves by accounting for station heights and negative 
antenna elevation angles, for several representative refractive indexes. 2 
Earlier extrapolations are based upon upper limits of the total bending 
associated with assumed earth-grazing rays. 3,4 

A refraction anomaly arising from a temperature inversion, storm, 
ducting, or other departure from an assumed representative radial ray- 
bending model precludes an absolutely confident evaluation of any given 
exposure at a given time. These phenomena are usually localized. How- 
ever, the intent of these computations is to protect the geostationary 
orbit against continuous interference arising simultaneously from a large 
number of terrestrial systems. 

Following a development of the refraction model, we derive the critical 
pointing azimuth corresponding to orbit intercept. The geocentered 
longitude displacement between the radio-relay site and the point where 
the refracted beam intercepts the orbit is next determined by spherical 
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geometry. Then the apparent slope* of the geometric orbit trace, as if 
viewed unbent by refraction from the station, is obtained corrected for 
aspect (also dependent upon antenna elevation angle and concomitant 
refraction). Using the refraction data converted to geometric elevation 
angles (as without ray bending or obstruction), the geometric orbit 
is adjusted to the apparent position and shape that would be observed 
at the radio site. This apparent orbit is hereafter termed the refracted 
orbit. Subsequent sections describe the determination of azimuth zones 
to be avoided for prescribed beam — orbit separations and the permis- 
sible transmitted power. Appendices are included to: provide means for 
estimating initial ray angles from commonly available radio-relay in- 
formation when actual antenna angles are unknown; solve by manual 
calculation a representative numerical example, giving the applicable 
equations for each step; and verify governing equations using a different 
analytical approach. 

II. DETERMINATION OF REFRACTION FOR POSITIVE AND NEGATIVE ANTENNA 
ELEVATION ANGLES 

2.1 Ray Tracing Equations 

In the following equations 6 is the initial angle* of a ray as it leaves 
the earth's surface and t is the total refractive bending corresponding 
to O . Figure 3 illustrates this relationship and shows the geometric 
director with elevation angle c to an intercept with the geostationary 
orbit at S. Also shown on an arc through S centered at P is the apparent 
position of the intercept S a and a horizontal reference, zero-elevation 
point A. The latter relationships are used in a subsequent section to 
describe a method for constructing the refracted orbit. 

Angle e is approximately, but in general not identical to, O — r*. 
However, this approximation is reasonable for rays between terrestrial 
antennas and geostationary satellites well beyond the earth's atmosphere 
when relatively small effects of parallax associated with the controlling 
portion of the atmosphere near the earth's surface are neglected. 6 

Figure 4 depicts a ray entering the earth's atmosphere and the result- 



* The first derivative with respect to pointing azimuth (azimuth-elevation plot 
of the orbit) is more completely defined in Section V. 

t 0. is generally used in ray-tracing equations to denote the initial ray elevation 
angle with respect to the local horizontal and is synonymous with or. used in 
subsequent sections to denote a radio-relay antenna beam elevation angle 
(namely, or in Figs. 6). # 

t Calculations using equation (5) show that the actual ray angle with respect to 
the geocenter differs from that resulting from the assumption e = 0„ — r by ap- 
proximately 15 minutes of arc for limiting conditions used herein. 
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Fig. 3 — Refracted ray. 

ing path of the ray due to the refractive gradient. Depending upon the 
angle of arrival of the ray when it enters the atmosphere, the refraction 
causes it to intercept the earth, graze the earth, or become tangent to a 
unique earth-centered sphere at some height above the surface. If the 
ray does not intercept the earth it continues out into space again, being 
subjected to approximately the same refraction in exit as it encountered 
upon entering. 

At any point on the ray trace, the angle the ray makes with the local 
horizontal at that point is denoted by 6. The angle between the tangents 
to the ray at any two points (a, b) is a measure of the refractive bending 
between them and is denoted by r(a, b). 

The following presentation of refraction is based largely upon equa- 
tions given in Refs. 2, 6, and 7 for zero or positive initial ray angles and 
as interpreted by the authors for application to negative angles. Appro- 
priate uses of the equations and their application to the problem are 
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Fig. 4 — Maximally refracted nongrazing ray. 
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explained, but the reader is directed to the references for complete 
derivations and limitations including sensitivities to the microwave 
frequencies involved. 

The radio refractive index n varies with atmospheric pressure, relative 
humidity, saturation vapor pressure and temperature. The lower limit 
for n is unity — no atmosphere, while the upper limit is determined by 
local climatic conditions. For the southeastern section of the United 
States, a typical range of n at mean sea level is 1.00025 to 1.0004. In 
equations involving refraction it is convenient to express refractivity 
as N where iV = (n — 1) X 10 6 or, for this case, N values are 250 and 
400 (N units). The average decay of N is approximately exponential 
with height and the difference in N between the surface and a height 
of 1 km above the surface is given by 2 

AN = -7.32 exp (0.005577iV.). (1) 

The subscript s denotes N at the surface. The decay constant with 
height is expressed by 2 

TV 

C = In — (2) 

' N. + AN K> 

N at any height h (kilometers above the radio site surface elevation) is 2 

N h = N.exp[-C e (h - h.)], (3) 

where h, is the surface elevation above mean sea level corresponding 
to AT, . 

N. is sensitive to local elevations and hence, charts of N, for moun- 
tainous regions are difficult to use because the N, contours are irregular 
and closely spaced. Therefore, obtaining an appropriate value of N. 
for use in equations (1), (2), and (3) for a particular site is often dif- 
ficult. However, charts are available giving N, reduced to mean sea 
level equivalents N„ which, in effect, reduce the height-dependent N, 
values to a common base. 8 Since charts of N are more easily interpreted 
and a single value of N usually applies over a large geographical area, 
they are used in this paper. 

N„ and N.(h) are related by N = N, exp (—h/7). Conversely, for a 
given value of N , N, for surface height h, above mean sea level is 
determined from the expression 

N. = N exp (-K/7), h. ^ 0. (4) 

The N, value obtained from equation (4) is used in equations (1), (2), 
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Fig. 5(a) — Radio refractivity conversion — N to N, (or Nr). 
Fig. 5(b) — 'Radio refractivity conversion — Nr to Nr- 

and (3). Figure 5a was generated using equation (4) and, for manual 
calculations, is convenient for determining N, . 

The conversion between N and N, above considers only the dry air 
effects of temperature and pressure related to the difference in elevation 
between mean sea level and the station surface height, whereas con- 
stants of the expression given for AN account for all terms contributing 
to a change in refractivity of the atmosphere with elevation above the 
surface height. 

A relationship exists between angle 0„ at the point of origin of the ray 
and at any other point on the ray trace which is expressed by 

(1 + 2V. X 10" 6 )(a + h.) cos e o 

= (1 + N h X 10" 6 )(a + h) cos d = C, (5) 

where C is a constant and a is the earth's radius at mean sea level (in 
kilometers). Thus, angle for any point on the ray trace is determined. 
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The incremental bending of the ray between closely spaced points 
on the trace is given by 

r(l,2) = -£'^COt0. (6) 

As suggested by Schulkin, the term 1/n is taken as unity with an error 
of less than 0.0001 in the computed refraction and, for an iterative 
solution, equation (6) is expressed as 7 

r An * 
r(l, 2) = - / cot BdAn. (7) 

Shulkin also shows 7 that equation (7) is approximated by (Awl — 
- An2)/0 m where An is n — 1 and m is (01 + 02)/2. Hence the in- 
cremental bending between two closely spaced points on the ray trace 
is expressed by 

Kl. 2) - ' m ' ~ *** 10 " 6 rad., S ,. S 10°, (8) 

0| -P 02 

where N x and N 2 , 6i and 2 are the N values and ray angles (in radians), 
respectively, at the closely spaced points. 

2.2 A Method for Calculating Refraction 

The following paragraphs describe a procedure for calculating re- 
fraction curves of the form in Figs. 6 for any values of N which is also 
applicable for the direct calculation of refraction corrections. 

2.2.1 Zero and Positive Initial Angles (+0 O ) 

First assume a value of N n , a station height h, , and an initial ray 
angle O . Equation (4) is used to determine N, for height h, and equa- 
tions (1), (2), and (3) to determine N for a height h, where h = h, + 
Ah*. Equations (5) and (8) then give the ray angle at the incremental 
height and the bending r in the first increment. For each successive 
increment of height, the previously solved-for values of N and become 
the initial values for equations (3) , (5) , and (8) . The values of r for each 
iteration are accumulated to give the total bending for h, and O . 
Repeating the above procedure with other values of O results in data 
points to be used in plotting the refraction curve for the assumed sta- 
tion height h. . A complete repetition of the above, beginning with 



* Incremental heights must be small in the lower atmosphere where n changes 
rapidly but may increase at the higher elevations. However, for the generation of 
Fig. 6, a constant increment of 0.25 km was used to a height of 90 km. 
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Fig. 6 — Refraction versus antenna elevation angle for (a) No = 250. and (b) 
V„ = 300. 

other station heights, results in a family of curves (namely, h,) for 
zero and positive initial ray angles. 

2.2.2 Negative Initial Angles ( — o ) 

The calculation of refraction for negative initial ray angles requires 
a modification of the technique used for positive angles. Note in Fig. 4 
that the bending of the ray from T3(/? 3 , — 3 ) to Tl is the same as that 
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Fig. 6 — Refraction versus antenna elevation angle for (c) N„ = 350, and (d) 
N. = 400. 
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of the ray from Tl (hi , 6 = 0) to T3. The latter is determined utilizing 
the equations for positive angles but in a slightly different manner. 

First, assume a grazing ray (h = 0, O = 0) and, using the iterative 
method previously described, compute the bending t and angle at 
each of the specific elevations h, desired for the chart, that is, in the case 
of Figs. 6:1,2, and 3 km. The value of t for each height is then added to 
the maximum t determined for the same height (when 0„ is zero as 
found in Section 2.2.1) to obtain the total bending for a ray with the 
initial angle, -0 O . In Fig. 4 this is T (T3, Tl) plus t(T1, S). Repeating 
the above but beginning with h = 1, 2, and so on, provides the data 
needed to extend the curves for positive initial angles into the negative 
range. 

Note in Figs. 6 that all refraction curves terminate on a dotted ex- 
tension of the zero-elevation ray representing a maximally refracted 
grazing ray. For a given height with an initial angle more negative than 
that represented by the point of termination, the ray is intercepted by 
the earth. 

2.2.3 Direct Machine Calculation 

For machine calculations it is desirable to compute the refractive 
bending directly without the use of refraction tables. For zero and 
positive initial ray angles, the calculation is straightforward as is des- 
cribed in Section 2.2.1. However, for a ray originating at a specific height 
with a negative initial angle O , it is first necessary to determine the 
height of the ray where is zero (h t at Tl in Fig. 4). A variation of 
equation (5) permits this determination. If the right side of equation 
(5) represents point Tl in Fig. 4 and the left the initial station, say T3 
at height h 3 , then 

(1 + 2V r3 X 10~ 6 )(a + h T3 ) cos 6 

= (1 + N T1 X 10 _6 )(a + h Tl ) = C. (9) 

The variables for the left side of equation (9) are all determined and 
evaluation yields the constant C. Then, values are assumed for h T1 
beginning with zero, N T1 determined, and the right side evaluated until 
the result is equal to C. After h Ti is determined, r(T3, Tl) and t(T1, S) 
are calculated as described in Section 2.2.2 and added numerically 
to give the total refractive bending. However, when h Tl is assumed to 
be for the initial iteration and the right side is larger numerically than 
C, the beam intercepts the earth and a solution with angle 0„ is not 
possible. In such cases it is often desirable to determine the initial graz- 
ing-ray angle to the mean sea level horizon, which is accomplished by 
incrementing 0„ upward until the equation is satisfied. 
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2.2.4 Angle from the Transmitting Site to the Radio Horizon 

The radio horizon, accounting for refraction and average local terrain, 
is determined by a variation of equation (9). Assuming that the height 
of a receiving station represents the average terrain height for some 
distance beyond, the initial angle a„ (= d ) of a ray which just grazes 
the terrain represents the angle to the radio horizon as seen from an 
elevated transmitting site. Letting the left side of equation (9) represent 
the transmitting location, the right side the receiving location (where 
is zero), and solving for a„ yields 



a H = cos 



(1 + N R X 10" ) (a + h 



_(1 + N T X 10-)(a + h 



r)J 



rad, h T ^ h R , a„ ^ 0, 

(10) 

where the subscripts R and T refer to the receiving and transmitting 
stations, respectively. N R is obtained by equation (4) and N T by equa- 
tions (1), (2), and (3), substituting N R and h R for N, and h, , and h T 
for h. (For manual calculations, N R and N T may be determined from 
Figs. 5. First enter Fig. 5(a) with h R , N , and read N R from the ordinate. 
Then enter Fig. 5(b) with (h T — h R ), N R and read N T from the ordinate.) 

2.2.5 Refractive Index Limits 

We now illustrate the importance of including the effects of refraction 
in orbital computations relating to terrestrial radio-relay systems. 
Table I (reflecting the use of equations and techniques discussed in 
subsequent sections) demonstrates parametrically the sensitivity to 
radio refractivity of the orbit-intercept pointing azimuth and the com- 
puted terrestrial transmitter power limitation. The 8 dB variation in 



Table I — Influence of Refraction 



Station Statistics 


Assumed Values 


Path azimuth 
Station latitude 
Station elevation 
Antenna elevation angle a 


103.5 degrees from true north 
55.0° north 
mean sea level 
degrees 


Parametric results 


Computed values 


Radio refractivity N, , N units 

Geometric elevation angle, e , degrees 
Critical azimuth (from north), degrees 
Maximum transmitter power, dBW 





250 


400 




102.6 

47 


-0.555 

101.76 

50.7 


-1.27 

100.7 
55 
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the power shown in the table indicates clearly that refraction must be 
included in calculations and that limits for refractivity should be care- 
fully considered. 

World charts of N in Ref. 8 indicate appropriate limits of N are 
250 and 400. At specific locations and for short time intervals, the index 
may not fall within these limits. However, localized conditions will 
not affect large numbers of stations at any given time and a wider 
range of N would unnecessarily broaden the restrictive zone for radio- 
relay systems. 

We suggest that the above limits be adopted for standardized cal- 
culations. For a specific case where an antenna pointing angle is close 
to the orbit and transmitter power limitations are restrictive, refractive 
limits applicable to that locale should be used. 

2.2.6 Adjustment of Computed Geometric Orbit Traces for Atmospheric 
Refraction 

For many solutions, particularly those involving graphical procedures, 
it is desirable to "elevate" a computed geometric orbit trace to its 
apparent (refracted) position and shape. Such an adjustment yields a 
presentation permitting a given, arbitrarily-shaped radio beam power 
profile to be related unrefracted and hence undistorted to the easily 
obtained configuration of the refracted orbit. Figures 7 are charts to 
enable this manipulation, produced from Figs. 6 by plotting (a — r a ) 
versus t„ . A method for using these charts is given in Section VI. 

III. DETERMINATION OF THE POINTING AZIMUTH TO ORBIT INTERCEPT 

Recall that the elevation angle (with respect to local horizontal) of 
the geometric director shown in Fig. 3 is denoted by e. Hence, station 
geometric elevation angle e„ may be replaced by a — r a0 where a is 
the initial antenna beam elevation angle and r a0 is the corresponding 
refraction correction. (A method for determining a„ is given in Appendix 
A.) 

Note that available information for established radio-relay routes 
in the United States giving antenna elevations, path distances, and 
antenna elevation angles is generally expressed in units of feet, statute 
miles, and degrees, respectively; it is necessary to convert these quanti- 
ties into kilometers and radians for use in many of the expressions which 
follow. 

Inspection of Fig. 8 shows that the azimuth displacement from the 
meridian through a station located at P to an intercept with the geo- 
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stationary orbit is identical to angle A of spherical triangle PES'. 
From laws for right spherical triangles 

. tan <p 
cos A = - — * , (11) 

tan B v 

where <p is the latitude at station P, and B is the arc equivalent of angle O 
of plane triangle OPS. Note that B is numerically equivalent to the 
maximum visible latitude for assumed radio-relay antenna beam eleva- 
tion angle a„ and concomitant total ray-bending angle r a0 . 

Angle B is determined from triangle OPS using the Law of Sines. 
This triangle is redrawn in Fig. 9, where 

sin 12 _ sin (t/2 + €„) 
a R 

(2 = sin" 1 (K~ % cose ), (12) 

where e„ is the station geometric elevation angle corresponding to a , 
a is the earth radius, and R is the orbit radius, R/a = K. From inspec- 
tion, = tt/2 — fi — e„ . Substituting for 12 yields 

jS = cos" 1 (A'" 1 cos «„) - e . (13) 

Substituting equation (13) for /3 in equation (11) results in 

A = cos" T . _ t .g?/ , -1 (14) 

Ltan [cos (K cos e ) — e„]J v 

IV. DETERMINATION OF RELATIVE LONGITUDE BETWEEN SITE AND ORBIT 
INTERCEPT 

The earth longitude displacement between radio-relay site P and 
suborbital intercept S' in Fig. S is side X of right spherical triangle 
PES'. From the Law of Sines 

sin X . 
— — t = sm 8, 
sin A 

from which: 

X = sin" 1 (sin A sin B), (15) 

where B and A are found by equations (13) and (14). Note that when 
A = tt/2, corresponding to <p = 0, the maximum visible longitude 
displacement is B. 
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Fig. 7 — Refraction versus geometric elevation angle for (a) No = 250, and 
(b) No = 300. 



V. GEOMETRIC ORBIT TRACE — CORRECTED FOR ANTENNA ELEVATION AND 
ATMOSPHERIC REFRACTION 

The geostationary orbit and earth's equator are coplanar; hence the 
orbit near the horizon normally appears to be tilted with respect to the 
local horizontal plane. Were an equatorial orbit sufficiently distant, 
as in celestial observations, the angle of tilt would equal the colatitude 
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Fig. 7 — Refraction versus geometric elevation angle for (c) iV„ = 350, and 
(d) N„ = 400. 
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of the observer's position. However, the maximum visibility circle in 
Fig. 8 is not a great circle, so the angle of tilt is less than the colatitude 
except where the latter circle crosses the equator. In the following para- 
graphs, an expression is developed for the tilt of the geostationary orbit 
as if viewed from a unique location on the visibility circle. This tilt is 
defined hereafter as the apparent slope of the orbit trace (or first deriva- 




Fig. 8 — Geometry relating site P and geostationary orbit intercept. 
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P 




R 
Fig. 9 — Geometry of intercept-determining triangle. 

tive of the orbit trace, in elevation, with respect to direction, in azi- 
muth).* 

The size of the visibility circle in Fig. 8 depends upon a and r aB ; 
for given site latitude <p, the angular displacement in earth longitude X 
between site P and point S where the refracted beam intercepts the 
geostationary orbit (longitude of suborbital point SO also depends upon 
a a and r ao . Hence, the slope 8 of the geometric orbit trace constructed 
in Fig. 11 also depends upon a and r ao . 

Note in Fig. 8 that the angle between orbit tangent t„ constructed at 
S and the local vertical plane at P through S is angle 4> of spherical 
triangle PES'. From the Law of Sines 



= sm-^- (16) 

\sin /3/ 

The complementary angle between orbit tangent t„ and the plane of a 
circle generated by radius CS (perpendicular to OP) is denoted by 
8' = tt/2 - 4>. Substituting equation (16) for <f> yields 

a*=cos-'(^Y (17) 

\sm /3/ 

Angle 5' is viewed in true magnitude from (or S'), but is seen from 
radio-relay site P as a smaller angle 8 when rotated through angle U 
as shown in Fig. 10. Note that tan 5' = y/x and tan 5 = (y cos 12) /x, 
from which 

8 = tan -1 (tan 8' cos Q,). (18) 

Combining equations (17) and (18) yields 

5 = tan" 1 {tan [cos -1 (sin <p/sm $)] cos G}. (19) 

* The orbit trace is envisioned as the locus of all pointing angles to_ the orbit, 
plotted on an azimuth-elevation chart aligned and calibrated according to the 
location of the observer. 
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An 




Fig. 10 — Development of geometric orbit slope as observed from site P. 



VI. CONSTRUCTION OF THE REFRACTED ORBIT 

Figure 11 presents the geometry of the problem viewed from a radio- 
relay station (similar, but not equivalent to the presentation given in 
Ref. 1). Reference to Figs. 1 and 3 may assist in interpretation of Fig. 11, 
wherein points A, S, and S ft correspond to similar points in Fig. 3 and 
intercept S corresponds to the left (easterly) intercept S in Fig. 1. Origin 
A is the beam-intercept direction from site P calculated from equation 
(14), accounting for atmospheric refraction. 

Since the orbit is tilted with respect to the local horizontal plane at 
site P, the elevation angle of the geometric director to the orbit con- 
tinuously varies as the orbit is scanned in azimuth. The refractive bend- 
ing of a ray is a function of the geometric angle, so that the position of 
the apparent, or refracted orbit, is above the geometric orbit and it 
exhibits a constantly changing slope with respect to the latter. 

The bent, refracted orbit is shown through point S a (apparent posi- 
tion of interception point S with refraction; also, the antenna elevation 
angle a at the azimuth origin). The straight line labeled "geometric 
orbit" is tangent at S (Fig. 8) to a radial projection of the geostationary 
orbit on a sphere of radius PS centered at P. The horizontal line shown 
through S also represents an arc, in edge view of a great circle through 
S, parallel at S to the local horizontal plane at P, on this same sphere. 
Hence, the angle 5 obtained from equation (19), corrected for refraction 
while retaining the concept of a constant site latitude <p, is also accurately 
represented by the apparent slope of the plane-figure geometric orbit 
trace at the azimuth origin. Figure 11 illustrates the construction of 
the corresponding refracted orbit trace. 

The equation of the linear geometric orbit trace through S is 



e' = —tan (5) A A + a — r ao , 



(20) 
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where e' is the elevation angle of the geometric orbit trace corresponding 
to an arbitrary displacement in azimuth A A from origin A. 

Figs. 7 are entered with values of d derived from equation (20) to 
obtain total refractive ray-bending angles r t > . The refracted orbit 
trace in Fig. 11 is then constructed by plotting points with coordinates 
(A + A -A, a!) and connecting these with a smooth curve, where 



= e' + TV 



(21) 



VII. AZIMUTH DISPLACEMENT FROM INTERCEPT TO KEEP THE BEAM CENTER 
AN ANGULAR SEPARATION V FROM (AND BELOW) THE MINIMALLY 
REFRACTED ORBIT (N MINIMUM) 

Figure 11 also illustrates a solution to keeping an angular separation 
v between the center of a circular beam and the geostationary orbit. 
The circle centered on S n and labeled "unrefracted beam" is a cross- 
section of a conical figure of revolution with apex at antenna site P 




"^\__j*EFRACTED 
BEAM 



Fig. 11 — Orbit geometry as observed from site P. 
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and 2v included angle (locus of all rays having angle v with respect to 
the beam center). If elevation angle a of the antenna is fixed by radio- 
relay path parameters, the orbit can be avoided only by an azimuth 
displacement of the beam. The required angular separation v between 
the orbit and the beam center results when the latter is moved in azi- 
muth away from intercept until the unrefracted cone is just tangent 
to the refracted orbit trace. Note that an identical relationship exists 
between the elongated "refracted beam" and the geometric orbit trace. 
However, since the former presentation is easier to develop and can 
readily accommodate unsymmetrical beam cross-sections, it is pre- 
ferred for following analyses. 

Figure 12 represents the analytical solution of an illustrative problem 
given in Appendix B and is helpful for visualizing subsequent pro- 
cedures. From inspection 

t«. a = (a ' ~ r - } ~ (a - = ° , (22) 

where a' is any point on the refracted orbit corresponding to arbitrary 
displacement M from an azimuth intercept. Note that in the range of 
interest M is negative with respect to azimuth A min , obtained from equa- 
tion (14) using e„ = a — r a0 . Within a few degrees of intercept, the 
geometric orbit is approximated by a line of constant slope 5. Recall that 
adjustment of the geometric orbit for refraction results in a refracted 
orbit trace having a constantly changing slope with azimuth displace- 
ment from the intercept. This displacement for a given elevation of the 
refracted orbit is found by solving equation (22) for M: 

M = («'-*-)- fa ~ O . (23) 

tan o 

Since the refracted orbit slope varies with elevation angle, no direct 
mathematical solution exists for determining that azimuth displacement 
providing an angular separation v between the beam center and the 
refracted orbit, measured in a direction normal to the latter. However, 
it is closely approximated by determining the slope of an orbit segment 
including the region of interest. Figure 12 suggests that the refracted 
orbit slope is everywhere less than the geometric slope. Hence, two 
judiciously chosen points on the refracted orbit having elevations 

al = a + V COS 8, a2 = a + V, 

always bracket the appropriate segment. Figures 6 and 12 also show 
that at these elevation angles the differential refraction is small. This 
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supports the assumption that the slope of the small orbit segment in 
this region is virtually constant. 

The refraction corresponding to elevation angles ct , al, and a2 are 
calculated or determined from the charts. Substituting a , t„„ , al, 
and T a i in equation (23) gives the azimuth displacement Ml corres- 
ponding to al. Similarly, M2 is determined by substituting a„ , t O0 , 
a2, and r a2 . Slope 51 between the two chosen points on the refracted 
orbit is given by 

tan 51 - («2 - al)/(M2 - Ml), 

resulting in 

51 = tan" 1 [v(l - cos 8)/ AM], (24) 

where AM = M2 - Ml. 

Figure 12 shows that azimuth displacement A A necessary to keep 
the beam center an angular separation v from the refracted orbit is 
S a c (equal to S a o plus oc).* By inspection, oc = v/sin 51, od = v/tan 51 
and S a o = M2 — od. Assembling these into an equation for A A yields 

AA miD = M2 - v/tan 51 + *>/sin 51. (25) 

7.1 Special Case 

If a ray with initial angle a intercepts the earth, a solution is obtained 
by calculating angle a n to the radio horizon as is described in Section II. 
Then, substituting e H for e in equation (14), the orbit intercept for a 
grazing ray is determined. When solving for the necessary azimuth 
displacement from this intercept, a H and r aH are substituted for a and 
r ao in equation (23). Note, however, that in determining al and a2 for 
substitution in equation (23), the use of ray angle a remains valid. 

VIII. AZIMUTH DISPLACEMENT FROM INTERCEPT TO KEEP THE BEAM 
CENTER AN ANGULAR SEPARATION V FROM (AND ABOVE) THE MAXI- 
MALLY REFRACTED ORBIT (iV MAXIMUM) 

The left side of Fig. 12 shows that the refracted orbit for the case 
of maximum refraction falls below the radio horizon with azimuth 
displacement from point A max obtained from equation (14) using e H = 
<xb — T aH . Since the earth intercepts all rays below the horizon, they 
cannot affect the orbit and it is only necessary to displace the beam 



* Notations such as S a c represent scalar distances between indicated points in 
Fig. 12. 
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center in azimuth sufficiently to maintain angular separation v from the 
horizon intercept. The required displacement is 

AA max = [v 2 - («. - a n f]K (26) 

In Fig. 12, a//(-0.25°) is the initial angle of a ray which grazes at a 
receiving station height of 0.4 km and originates at an assumed trans- 
mitting height of 0.5 km. 

As demonstrated in Appendix B, little increase in the critical zone 
results if it is assumed that the angle to the radio horizon is equal to 
the initial ray angle. Therefore, for manual calculations involving small 
positive values of a , these angles are assumed identical so that AA 
is simply an angular separation v from the orbit intercept determined 
for a . For values of a e more negative than that for a grazing ray 
(determined from Figs. 6 or by calculation), it is necessary to determine 
a H and the corresponding orbit intercept using equations (10) and 
(14). Then equation (26) gives the necessary azimuth displacement. 

IX. DETERMINATION OF THE CRITICAL ZONES 

9.1 Critical Zones Defined 

The critical zones to be avoided at radio-relay transmitting sites to 
protect the geostationary orbit are defined: 
ZcHt = A max + AA max to A miD - AA min 

(degrees from South), (27) 

where values for A and AA for maximum and minimum refraction are 
obtained as in Sections III, VII, and VIII. These zones are converted 
to azimuth zones with respect to true north by subtracting the bound- 
aries from ISO degrees for the easterly zone and adding them to 180 
degrees for the westerly zone. 

Calculations for stations in southern latitudes are identical except 
that they are referenced to north rather than south. The easterly azi- 
muth zone with respect to true north is obtained directly from equation 
(27), while the zone boundaries are subtracted from 360 degrees for the 
westerly zone. 

9.2 Special Case 

At latitudes exceeding <p = cos -1 (K~ l cos \^) — ^ it is impossible 
for the beam's center ray to be below the orbit with angular separation 
v.* Hence, for such extreme northern latitudes, a single critical zone 

* Derivation of this equation is given in Appendix C. 
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spans south with both easterly and westerly boundaries denned by 
Z CTit = A max + AA max (degrees from south). (28) 

X. DETERMINATION OF MAXIMUM PERMISSIBLE RADIATED POWER 

As mentioned in Section I, international agreements exist for maxi- 
mum radiated powers. 9 For 6-GHz radio-relay transmitters whose 
antennas point within v = 2° of the geostationary orbit, the power 
limitation for separations less than 0.5° is 47 dBW relative to the iso- 
tropic case (EIRP), increasing 8 dB per degree to a maximum of 55 dBW 
(occurring at 1.5°). This limitation refers specifically to the center of 
the major lobe. For this case, only the relationship between the refracted 
center ray of the beam and the geometric orbit is considered. 

If an existing or proposed path direction is between the critical values 
computed as in Section III for maximum and minimum refraction, the 
refracted center ray is likely to intercept the orbit for appreciable periods 
of time. For such cases the maximum power is 47 dBW. 

If the path direction for a system in the northern hemisphere is within 
the critical zone but nearer to south (smaller) than A mia , the actual 
separation p illustrated in Fig. 12 is 

p = (A mia - A p ) sin 5, (29) 

where A v is the path direction measured from true south. 

Conversely, if A P exceeds A max , the separation of the beam center 
from intercept of the geometric orbit and the refracted horizon (Fig. 12) 
is 

P = [(A v - A max ) 2 + | e* - 6 | 2 )*, (30) 

where 

*h = ctu — r aH , e = a ~ T ao . 

Note that e in equation (30) represents maximum atmospheric refrac- 
tion. 

If a ray having initial angle a intercepts the earth, e„ for use in equa- 
tion (30) is indeterminate. For such special cases, a conservative ap- 
proximation for the angular separation is p = A p — A mBX . 

The maximum permissible effective radiated power P t dBW (EIRP) 
is given for separation p according to the following criteria: 

P ^ 0.5°, P, = 47; 
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0.5° < p ^ 1.5°, P, = 8(p - 0.5) + 47; (31) 

p > 1.5°, P t = 55. 

XI. CONCLUSIONS 

A direct analytical method involving few approximations and assump- 
tions can be used by system planners for calculating refraction-cor- 
rected ranges of pointing azimuth for microwave radio-relay antennas 
within which significant interference with geostationary communication 
satellites can be expected. Required angular separations between the 
refracted beam and the geostationary orbit are translated into required 
azimuth displacements of a radio-relay antenna from that calculated 
for orbit intercept; conversely, for cases where exposure is unavoidable, 
means for determining the maximum transmitted powers permitted 
by international agreement are presented. 

Since all analytical expressions including refraction corrections are 
readily amenable to machine calculation, both speed and improved 
accuracy in estimating the pointing azimuths are possible. The sug- 
gested refractive index limits are believed to be representative for the 
large majority of exposures and useful for a standardized approach to 
the problem. For more general applications, where refractive index 
variations are known to be different, one may use the same principles 
to generate his own applicable correction curves. 
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APPENDIX A 

Estimation of Antenna Elevation Angles 

The initial beam elevation angle for a radio-relay antenna is deter- 
mined by the geometry of transmitting and receiving locations, the path 
length, and refraction. The final alignment based upon transmission 
measurements, if recorded, is preferred for these calculations. The 
antenna elevation angle for proposed radio-relay paths can be estimated 
using the method given below with sufficient accuracy. 
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Figure 13 depicts radio-relay transmitter T and receiver R at eleva- 
tions OT and OR above geocenter 0, assuming a spherical earth of radius 
ka. Coefficient k is the ratio of apparent earth radius to true earth radius 
and accounts for refraction in the lower atmosphere. 10 The path length 
is represented by arc D (great circle length at mean sea level). The trans- 
mitting and receiving antenna heights with respect to mean sea level 
are denoted by h T and h R , respectively. 

Inspecting Fig. 13, <£ = D/ka radians, C = 2(ka + h T ) sin (<f>/2), 
and a = E — <f>/2 radians. From triangle TT'R and laws for plane 
triangles 



t&nE = 



(ha - M sin fr/ 2 ~ ft/ 2 ) 



C + (h R - h T ) cos (tt/2 
It can also be shown that 



ft/2) 



(32) 



a„ = tan 



■[; 



;]- 



D/2ka radians, (33) 



-tan (D/2ka) X (2ka + h R + h T ). 

where h R , h T , D, and a are expressed in the same units. 

Reference 2 provides a formula and a table relating k and N a . For 
most calculations, a value of k = 4 results in sufficient accuracy 10 (the 



TT/2 + <t>jz 




ka+hi 



Fig. 13 — Geometry of antenna elevation angle. 
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geometries of radio-relay systems within limits of normal site elevations 
require antenna elevation angles which are relatively insensitive to 
values chosen for k). 

APPENDIX B 

Illustrative Calculation 

Problem Input for a Particular Radio-Relay Site: 

Latitude <p — 38°N, 

Transmitter Height h T — 0.5 km, 
Receiver Height h R — 0.4 km, 
Path Length D — 28 km, 

Path Azimuth A v — 97.75° with respect to true north, equivalent 

to 82.25° from south towards east, 
N Limits — 250 and 400 N units. 

B.l Antenna Elevation Angle 

From equation (33) and using k = | 

f" 0.4 - 0.5 1 

_tan [(28/(2.66 X 6373)](2.66 X 6373 + 0.4 + 0.5) J 

28 



a n = tan" 



2.66 X 6373 
= -0.00523 rad, which converts to —0.3°. 

B.2 Geometric Elevation Angle jor N„ = 250 
From Fig. 6a 

e o = a „ - r Q „ = -0.3 - 0.59 = -0.89°. 

B.3 Azimuth Intercept for N„ = 250 
From equations (13) and (14) 

= cos" 1 [0.1509 cos (0.89°)] + 0.89 = 82.2°, 
-4 min = cos" 1 [tan (38°)/tan (82.2°)] = 83.86° from south. 

b.4 Orbit Slope jor N = 250 
From equations (12) and (19) 
Q = sin -1 [0.1509 cos (0.89°)] = 8.68°, 
8 = tan" 1 [tan [cos" 1 [sin (38°)/sin (82.2°)]} cos (8.68°)] = 51.26°. 
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B.5 Azimuth Displacement for v = 2° below the Orbit 
From equations (23), (24), (25), and Fig. 6(a) 

a = -0.3°, r ao = 0.59°, 

a\ = -0.3 + 2 cos (51.26°) = 0.95°, r al = 0.35°, 
a 2 = -0.3 + 2 = 1.7°, r a2 = 0.28°, 

Ml = [(0.95 - 0.35) - (-0.3 - 0.59)] + tan (51.26°) = 1.19°, 
M2 = [(1.7 - 0.28) - (-0.3 - 0.59)] + tan (51.26°) = 1.85°, 
51 = tan" 1 [2[1 - cos (51.26°)]/0.66| = 48.62°, 

AA min = 1.85 - 2/tan (48.62°) + 2/sin (48.62°) 
= 2.76° toward south from intercept. 

B.6 Horizon Intercept for N = 400 
From equations (1), (2), (3), (4), (10), (12), (13), (14), and Fig. 6(d) 

N R = 400 exp (-0.4/7) = 377.78 or from Fig. 5(a), 

AN = -7.32 cxp (0.005577 X 377.78) = -60.2, 

C. - In [377.78/(377.78 - 60.2)] = 0.17, 

N T = 377.78 exp [-0.17(0.5 - 0.4)] = 371.28 or from Fig. 5(b), 

., f l. 000377 6373.4 "! 
a « = - C0S LrOOOSTl X 6373^J = ~ ' 25 ' 

r aH = 1.22°, 
e„ = -0.25 - 1.22 = -1.47°, 
= cos -1 [0.1509 cos (1.47°)] + 1.47 = 82.79°, 
A max = cos -1 [tan (38°)/tan (82.79°)] = 84.33° from south. 
b.7 Azimuth Displacement for v = 2° from Horizon Intercept 
From equation (26) 

AA max = [(2) a - (-0.3 + 0.25) 2 ]* 

= 1 .99° toward north from intercept. 

b.8 Critical Zone 
From equation (27) 
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Z eTit = 84.33 + 1.99 to 83.86 - 2.76 
= 86.3° to 81.1° from south. 

B.9 Azimuthal Zones 
True-north azimuths: 

93.7° to 98.9° easterly, and 
266.3° to 261.1° westerly. 
The path azimuth of 97.75° falls within the easterly azimuthal zone. 
B.io Relative Longitude to Suborbit Intercepts 
From equation (15) 

X min = sin" 1 [sin (83.86°) sin (82.2°)] = 80.04°, 
X ma x = sin -1 [sin (84.33°) sin (82.79°)] = 80.82°. 

B.n Maximum Permissible Radiated Power 

Comparing the path direction of 82.25° from south with the critical 
zone found in Section B.8 and with A mU and A max found in Sections B.3 
and B.6 reveals that equations (29) and (31) are appropriate for calculat- 
ing the angular separation and permissible power: 

P = (83.86 - 82.25) sin (51.26°) = 1.23°, 
P, « 8(1.23 - 0.5) + 47 = 52.8 dBW at 6 GHz. 

B.ii.i Alternate Power Calculation 

An alternate calculation is indicated when the path direction is further 
from south than A max . Assume that A„ is, instead, 85.25° from south 
(left side of Fig. 12). The appropriate equations are now (30) and (31). 
A value of e„ for maximum refraction is required for equation (30). 

From Section B.5, Fig. 6(d), equations (30) and (31) 

a = -0.3°, 

T a . = 1-3°, 

€„ = -1.6°, 

p = [(85.25 - S4.33) 2 + (-1.47 + 1.6) 2 ]* = 0.93°, 
P t = 8(0.93 - 0.5) + 47 = 50.4 dBW at 6 GHz. 
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APPENDIX C 

Alternate Derivations of Basic Equations and Critical Latitudes 

All following relationships are expressed in terms of latitude ip of radio- 
relay site P v shown in Fig. 14 (a constant), and the maximum latitude 
for which the refracted geostationary orbit is visible to the antenna 
(a constant elevation angle). 
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Fig. 14 — Anglp <f>i between orbit tangent and local vertical plane 0-S-P v . 
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c.i Latitude Limit of Visibility for the Geostationary Orbit 

Figure 14 shows that the maximum latitude visible to the geosta- 
tionary orbit is 



= cos ' [K cos (a — r ao )] — a — 



T, 



(34) 



= /3, 
where 



a = known elevation angle of the radio-relay antenna, 
r ao = corresponding total ray-bending angle due to refraction in- 
ferred from CRPL Exponential Reference Atmosphere* or from Figs. 6 
of text, and 

K = ratio of orbit radius R to assumed earth radius a. 

Hence, equation (34) is a restatement of equation (13). 

C.2 Pointing Azimuth to Orbit Intercept 

The larger of angles formed by intersection of site latitude circle <p 
and the visibility circle corrected for antenna elevation and refraction 
shown in Fig. 14 is a direct measure of the pointing azimuth to be avoided 
for a station located at that intersection (P v ). Hence, angle A N between 
tangents T p and T v is the supplement of pointing angle A given by 
equation (14) referred to true south. From the geometry of Fig. 14 

A N = cos -1 S — sin tan -1 (7-7: r^ 2 .» , ) 

[ L M (1 — sm ip — cos py |/_ 

-h-'( l(1 - sin :;r°^ M )l 

tt/2 = I A N I < r. (35) 

Equation (35) yields two pointing azimuths which are of interest; 
one in the second quadrant referred to true north, corresponding to a 
westerly direction for stations in the northern hemisphere — and one 
in the third quadrant, or easterly direction. Reversed directions result 
for stations in the southern latitudes. 

Figures 15 illustrate changes of variables which simplify the demon- 
stration of equivalence between equation (35) and equation (14). Since 
cos 2 = 1 — sin 2 /3 and cos 2 tp = 1 — sin 2 <p, 



A N = cos \ — s 



. r. -1 / sing \1 

sm tan [■ , , ■ ■> Q .—5 — vrij 

L \| (sm - sm <py \/J 
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COS 



tan" 1 f 



I (cos 2 g - cos 2 ft)* I 
cos /3 



Substituting n and v (Figs. 15) and reducing the result yields 

Aft = cos -1 (—sin n cos v). 
Now, substituting for n and p provides 

A N = cos -1 [— (sin (p/sin /3)(cos /3/cos <^)] 

tan yO 
tan /3 



= cos l — 



which is exactly the supplement of the angle obtained from equation (11) 
in the text. 

c.3 Longitude Displacement of Orbit Intercept 

The earth-longitude displacement between radio-relay site P v and 
suborbital intercept S' in Fig. 14 is also inferred from equation (35): 



rl /, -2 2 
(1 — sin to — cos 
J_i £ 



X = tan 



cos /3 



ts\ 



^ I X I < ir/2, 



= sin ' (sin A sin 0), (36) 

from which first and fourth quadrant longitude adjustments are ref erred 
to the suborbital longitude. The equivalence with equation (15) is 
demonstrated using techniques indicated above and identifying angle /3 
uniquely with the maximum latitude for visibility ip^ (Section C.l). 

c.4 Geometric Orbit Trace 

Since the size of the visibility circle in Fig. 14 depends upon a and 
T ao , the earth-longitude displacement X between P v and S (same 
longitude as suborbital point S') also depends upon a and r ao . Because 




SIN 9 




(cos^v-ccs 2 /?) 2 ! 



cos/3 



(SIN 2 /?-SIN 2 ?>)2| 

(a) (b) 

Fig. 15 — Change of variables: (a) angle /* and (b) angle v. 
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the resulting viewing aspect depends upon these parameters, the slope 
of the geometrical orbit trace constructed in Fig. 11 also depends upon 
a and r ao as well as <p. 

Note that the angle between orbit tangent t constructed at S and the 
local vertical plane at P v through S is denoted by 4> t . This angle is 
observed undistorted at point S', but appears to be a slightly enlarged 
angle <f>' t for an observer at P v (as if point S were visible without ray 
bending) : 

0{ = tan -1 (tan 0,/cos $2). (37) 

Hence, complementary angle 5 between orbit tangent t and a line 
through S perpendicular to the geometrical line-of-sight SP V and parallel 
to the horizontal plane at P v is also the slope of the corrected geometric 
orbit trace shown in Fig. 11, 

8 = cot" 1 (tan </>,/cos £2) 

tan [sin" 1 (sin <p/sin /3)] 



= cot" 



= tan" 



cos {tan -1 [sin p/(K - cos 0)]\ 
{tan [cos -1 (sin cp/sm /3)] cos £2} , 



(38) 

for which equivalence to equation (19) is shown using the techniques 
incorporated in Sections C.2 and C.3. 

C.5 Maximum Latitude Permitting Angular Separation v Below the Orbit 

A maximum latitude exists for each antenna elevation angle allowing 
the beam-center ray to be below the refracted orbit with prescribed 
angular separation v. Figure 16 illustrates the determination of this 
critical latitude. S p is a point on the geostationary orbit having zero 
relative longitude with respect to station site P v . A ray emanating from 
the antenna with initial vertical angle (tt/2 + a + v) must just inter- 
cept the orbit. Accounting for refraction, the angle between the geo- 



«A=a +y-T (ao + y .) a -T a o p v. 




R I ~0 

Kip;. 16 — Geometry determining the critical latitude. 
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metric director of the elevated ray and the local horizontal is ^ = 
[a + v — T (ao+V) ]. Site latitude <p and M> are related by triangle OP v S p . 
From the Law of Sines 

sin (tt/2 + fr) sin fl 
R a 

Letting R/a = K, 

fi = sin" 1 (K' 1 cos ¥). 

Since fi + <? + ir/2 + ^ = t, 

<p = cos -1 (K _1 cos ¥) - *. (39) 
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